Photon emission by an ultra-relativistic particle 
channeling in a periodically bent crystal J 

Andrei V Korolf^ d, Andrey V Solov'yov| |, and Walter Greiner % 

fDepartment of Physics, St. Petersburg State Maritime Technical University, Lcninskii 
prospect 101, St. Petersburg 198262, Russia 

lA.F.IofFe Physical- Technical Institute of the Academy of Sciences of Russia, 
Polytechnicheskaya 26, St. Petersburg 194021, Russia 

^Institut fiir Theoretische Physik der Johann Wolfgang Goethe- Universitat, 60054 
Frankfurt am Main, Germany 

Abstract. This paper is devoted to a detailed analysis of the new type of the 
undulator radiation generated by an ultra-relativistic charged particle channeling along 
a crystal plane, which is periodically bent by a transverse acoustic wave, as well as to 
the conditions limiting the observation of this phenomenon. This mechanism makes 
feasible the generation of electromagnetic radiation, both spontaneous and stimulated, 
emitted in a wide range of the photon energies, from X- up to 7-rays. 



PACS numbers: 41.60 
1. Introduction 

This paper is devoted to a detailed analysis of the new type of the undulator radiation 
(acoustically induced radiation - AIR) generated by an ultra-relativistic charged particle 
channeling along a crystal plane, which is periodically bent by a transverse acoustic wave 
(AW), as well as to the conditions limiting the observation of this phenomenon. This 
mechanism, suggested recently in , makes feasible the generation of electromagnetic 
radiation, both spontaneous and stimulated, emitted in a wide range of the photon 
energies, from X- up to 7-rays. 

The mechanism of the AIR generation is illustrated in figure Under the action of 
a transverse acoustic wave propagating along the z-direction, which defines the center 
line of an initially straight channel (not plotted in the figure) the channel becomes 
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periodically bent. Provided certain conditions are fulfilled (see below in this section and 
in section |^), the beam of positrons, which enters the crystal at a small incident angle 
with respect to the curved crystallographic plane, will penetrate through the crystal 
following the bendings of its channel. It results in the transverse oscillations of the 
beam particles while travelling along the z axis. These oscillations become an effective 
source of spontaneous radiation of undulator type due to the constructive interference 
of the photons emitted from similar parts of the trajectory. As we demonstrate below, 
the number of oscillations can vary in a wide range from a few up to a few thousands 
per cm depending on the the beam energy, the AW amplitude, a, and wavelength. A, 
the type of the crystal and the crystallographic plane. In addition to the spontaneous 
photon emission by the undulator, the scheme presented in figure leads to a possibility 
to generate stimulated emission. This is due to the fact, that photons, emitted at the 
points of the maximum curvature of the trajectory, travel almost parallel to the beam 
and, thus, stimulate the photon generation in the vicinity of all successive maxima and 
minima of the trajectory. 
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Figure 1. Schematic representation of the initiaUy hnear planar channel bent by the 
transverse acoustic wave. The notations are: d is the channel width, a, A are the AW 
amplitude and wavelength, respectively. 

Some aspects of the electromagnetic radiation by a beam of charged particles 
channeling in a crystal undergoing the action of an ultrasonic wave have been discussed 
in the literature We shall mention the main results obtained in the cited papers 

and outline the features, both of formal and of principal nature, which distinguish the 
AIR radiation described above from the effects considered in 0-0 at the end of this 
section. 

Prior to that let us briefly discuss the phenomena closely related to the mechanism 
of generating the acoustically induced radiation. These are the channeling effect, the 
channeling radiation and the undulator radiation. 

The basic effect of the channeling process in a straight crystal is that a charged 
particle can penetrate an anomalously large distance inside a crystal if travelling nearly 
parallel to the crystallographic plane or axis and experiencing the collective action of 
the electric field of the lattice ions. Positively charged particles are steered into the 
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interatomic region (the planar channeling) while negatively charged projectiles move in 
vicinity of the ionic chains (the axial channeling). 

Channeling was discovered in the early 1960s by computer simulations of ion motion 
in crystals |T0| , |Tl]]. Large penetration lengths were obtained for ions incident along 
crystallographic directions of low indexes. It is interesting to note that such a guided 
motion had already been predicted |T^ by Stark in 1912 [0. A comprehensive 
theoretical study |T^, |TB[ introduced the important continuum approximation for the 
interaction potentials between energetic charged projectiles and lattice atoms arranged 
in atomic strings and planes. These concepts were subsequently widely used in the 
interpretation of channeling experiments (see e.g. the review article ll^). During recent 
years particular attention has been paid to studying of the channeling phenomenon at 
high energies of the particles |]16|, |17[ 0). 

The idea of bending the high-energy beams of charged particles, proposed in ||19|| , 
has become of great current interest 



16| , pTSf , HH because of its practical 



application for the manipulation of the beams. Indeed, the high-energy beams can be 
bent by crystals as efficiently as by means of external electric or magnetic fields. The 
size of such crystals is typically hundreds times smaller than the size of the equivalent 
magnet. In recent experiments with 450 GeV protons p3| the efficiency of the particle 
beam deflection was reported on the level of 60%. 

Another possible application of bent crystal was suggested in ||2^ , where the question 
whether the channeling effect can also be used to focus beams, especially heavy-ion 
beams, was discussed. To this end a crystal is needed in which the crystal axes are no 
longer parallel, but are slanted more and more the farther away they are from the axis 
of the beam. Then the bending angle of the particles far away from the beam axis would 
be largest and a general focusing effect will result. Such a crystal can in principle be 
produced by varying the germanium to silicon ratio in a mixed crystal |^ . 



The criterion for a stable channeling of an ultra-relativistic particle in a bent crystal 
was formulated in |19| and has clear physical meaning: the maximum centrifugal force 
due to the channel bending must be less than that of the interplanar fleld. Provided 
this condition is fulflUed, the beam of (positively charged) channeling particles at each 
instant moves inside the channel, especially parallel to the bent crystal midplane as it 
does in the case of linear channeling. 

Two realistic ways of "preparation" of a periodically bent channel may be discussed. 
It is feasible, by means of mordern technology [25|, to grow the crystal with its channels 
been statically bent according to a particular pattern. The channeling phenomenon 
has been also considered in the periodic structures such as superlattices 0, 0. THe 
possibility fl], |2[, p - [Q, of a dynamically bent crystal, arises if one considers 
propagation of an acoustic wave along some particular direction in a crystal as shown 
in flgure |I|. Under the action of the AW the crystal channels, being linear initially. 
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will be periodically bent. In both cases the passage of an ultra-relativistic particle 
along the bent channel gives rise to the AIR phenomenon due to the curvature of the 
trajectory, provided that the projectile is trapped into the channel. The latter condition 
is subject to the general criterion for the channeling process in a bent crystal ||19| (see 
also 1^, and can be fulfilled by a proper choice of projectile energy and maximal 

curvature of the channel, as described below in the paper. 

The advantage of the static channel is that its parameters are fixed and, thus, the 
projectile moves along the fixed trajectory as well. To calculate the characteristics of 
the emitted radiation one needs to know only the number of the periods and the local 
cuvature radius. The disadvantage is that when fixing the number of periods then the 
parameters of the system can be varied only by changing the energy of the particle. 
This makes the photon generation less tunable. 

The motion of a particle in the field of the channel bent by transverse AW becomes 
periodic and, therefore, the spectrum of radiation emitted by a projectile acquires 
essential features of the undulator radiation. 

It is natural that the AIR emission is accompanied by the ordinary channeling 
radiation p8|, |2^. This specific type of electromagnetic radiation arises due to the 



transverse motion of the channeling particle inside the channel under the action of the 
interplanar field. The phenomena of the channeling radiation of a charged projectile in 
a linear crystal (see eg. [^, ^) as well as in a "simple" (i.e. non-periodically) bent 



channel |3^, are known, although in the latter case the theoretical and experimental 
data are scarce, at least up to now. 

We demonstrate in our paper that the AIR and the ordinary channeling radiation 
can be separated provided the condition a 3> (i is fulfilled (see figure p. Then the 
frequency of the particle oscillations inside the channel is much higher than the frequency 
of the transverse oscillations caused by the AW. Therefore, these two motions are well 
separated and the AIR mechanism can be treated independently from the ordinary 
channeling radiation. We shall discuss this important issue in more detail in section 
Here let us only note that a similar situation occurs in a one-arc bent crystal, where the 
channeling charged particle generates additional synchrotron type radiation due to the 
curvature of the channel [^, 33]. This component of the total radiation intensity leads 
to the undulator effect in the channel periodically bent by AW. We will show that the 
intensity of AIR can be made larger than that of the channeling radiation. 

As was pointed out in [jl], 0, the system "ultra-relativistic charged particle + 
periodically bent crystal channel" represents by itself a new type of undulator, 
and, consequently, serves as a new source of undulator radiation of high intensity, 
monochromaticity and with a particular pattern of the angular-frequency distribution. 
The electromagnetic radiation in this undulator arises mainly due to the bending of 
the particle trajectory which follows the shape of the channel. The parameters of this 
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undulator, as well as the characteristics of the electromagnetic radiation, depend on the 
type of the crystal and on the crystallographic plane, on the type of projectile and its 
energy; they also depend on the shape of the bent channel, and, thus, can be varied 
significantly by varying the enlisted characteristics. 

The theory and also the various practical implementations of the undulator 
radiation, e.g. the radiation emitted by a charge moving in spatially periodic static 
magnetic fields (a magnetic undulator), in a laser field (a laser-based undulator), etc. 
have a long history [03, IH^ and are well elaborated PUI IHT], EO, P7[ . It is also known 



that a positively charged ultra-relativistic particle, undergoing a planar channeling in a 
linear crystal, radiates electromagnetic waves whose spectral and angular distributions 
are those of the undulator (a natural undulator) due to the transverse oscillations caused 
by the action of the repulsive interplanar potential. Also the characteristics of the 
radiation emitted in a natural undulator are very sensitive to the distribution of the 
particles in the beam over the transverse energy and the incident angle. 

The important feature of the dynamically bent crystal by means of an AW is that it 
allows to consider an undulator with the parameters and p (here is the number 
of periods in the undulator and p = 2717 a/ A is its parameter (see figure |l|), and 7 
is the relativistic factor) varying over a wide range, which is determined not only by 
the projectile's energy but also by the AW frequency and amplitude. The latter two 
quantities can easily be tuned resulting in the possibility of significantly varying the 
intensity and shape of the angular distribution of the radiation. It is important to 
note that the parameters of the acoustically based undulator are inaccessible in the 
undulators based on the motion of charged particles in periodic magnetic fields and also 
in the field of laser radiation |^, ^|. 

In the suggested scheme AIR is generated by the relativistic charged particles, with 
relativistic factors 7 = e/mc^ ^ 1 (c is the velocity of light, m is the mass of the particle 
and e is its energy). The large range of 7 available in modern colliders for various 
charged particles, both light and heavy, together with the wide range of frequencies 
and the amplitudes possible for AW in crystals allow to generate AIR photons with the 
energies up into the GeV region [||]. 

As demonstrated in 0, the specific pattern of the undulator radiation combined 
with the AIR mechanism allows to discuss the possibility to create a powerful source 
of stimulated monochromatic radiation of high-energy photons. We analyse this 
important possibility in more detail in the present article. Particular attention is 
paid to the investigation of the infiuence of various physical processes, such as the 
photon attenuation and the dechanneling effect of positrons in crystals, both leading to 
a reduction in the amplification of the stimulated photon emission. 

Let us briefiy outline here the assumptions, which are adopted in this paper and 
formulate the main conditions which allow to consider the AW bent crystal as an 
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undulator. 

Although the general treatment of both the undulator radiation (section 0) and the 
channeling process in the AW bent crystal (section ^) can be applied to axial and a 
planar channeling, we consider the latter case because it is known (see e.g. that 
for bent crystals negatively charged projectiles are steered along crystallographic axes 
much less efficiently than positively charged ones along the crystal planes. In fact, 
experimentally the axial channeling of negative particles has not been observed so far. 

Two important conditions which we assume to be fulfilled are the following: 

a) The bending of the channel becomes significant if the AW amplitude is noticeably 
larger than the interplanar spacing, d (see figure p. Thus, the strong inequality 

a > (i (1) 

is implied throughout the paper. Typically, the spacing between the planes, which 
are characterized by the low values of the Miller indices (such as the (100), (110) 
and (111) planes), lies within the range 0.6. . .2.5 A(see e.g. [^] and table |l| in 



section H) below. Therefore, (|l]) is fulfilled for the AW amplitudes a > 10 A. 

b) The channeling process in a bent crystal takes place if the maximal centrifugal force 
in the channel, m7f^/i?mm (where -Rmin is a minimum curvature radius of the bent 
channel) is less than the maximal force due to the interplanar field |l^, |^, ^ : 



m7T;Vi?min < gf/Lx- (2) 

Here q is the charge of the projectile and the quantity t/^ax stands for the maximum 
gradient of the interplanar field. 

Provided (^ is fulfilled, the projectile, incident at small angle with respect to the 
mid-plane, will be trapped in the channel (see figure [^). The condition (||) allows 
to disregard the oscillations of the particle inside the channel due to the interplanar 
force q U' . Therefore, the particle trajectory will be primarily harmonic defined by the 
equation 

?/(z) = asin (^27r^^ , ^ = [0 . . . L] (3) 

Here L stands for the thickness of the crystal. The minimum curvature radius of this 
trajectory is equal to -Rmin = (A/27r)^/a. Thus, the decrease in -Rmin and, consequently, 
the increase in the maximum acceleration of the particle in the channel is achieved 
by decreasing A and increasing a. As a result, photon emission due to the projectile's 
acceleration in the bent channel may be significantly enhanced. This radiation is emitted 
coherently from similar parts of the trajectory, and may dominate |1[] over the channeling 
radiation. 
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An adequate approach to the problem of the radiation emission by an ultra- 
relativistic particle moving in an external field was developed by Baier and Katkov in 
the late 1960s and was called by the authors the "operator quasi-classical method" . 
The details of that formalism can be found in |31|, ^ and will not be reproduced here. 

The advantage of this method is that it allows to use the classical trajectory for the 
particle in an external field and, simultaneously, it takes into account the effect of the 
radiative recoil. 

The classical description of the particle motion is valid provided the characteristic 
energy of the projectile in an external field, hujo-, is much less than its total energy, 
e = m^c^. The relation huo/e oc 7"^ ^ 1 is fully applicable in the case of an ultra- 
relativistic projectile. Typical values of Tiojo are in the order of magnitude equal to the 
continuous potential-well depth, vary from several eV for the crystals of light materials 
(e.q. the LiH crystal |Q) up to 10^ eV for the crystals of heavy materials (e.g. the W 
crystal |2|). 

The role of the radiative recoil, i.e. the change of the projectile energy due to 
photon emission, is governed by the ratio %uj/e. In the limit Tluj/e a. purely classical 
description of the radiative process can be used. For Tiuj/e < 1 the quantum 
corrections due to the emission of the photon must be taken into account. 

The quasi-classical approach neglects the huo/e terms, but it explicitly takes into 
account the quantum corrections due to the radiative recoil in the whole range of the 
emitted photon energies, except for the extreme high energy tail of the spectrum. 

Using this method the spectra of photons and electron-positron pairs in linear 
crystals were successfully described It was also applied to the problem of 

synchrotron-type radiation emitted by an ultra-relativistic projectile channeling in a 



non-periodically bent crystal ^ 



We use this general formalism of Baier and Katkov in our paper to treat the AIR 
phenomenon. 

To conclude the introductory section let us discuss the results of the papers [Q- 
where the electromagnetic radiation emitted by a charge passing through a crystal 
undergoing the action of an ultrasonic wave was considered. 

As far as we are aware of, the first study of the infiuence of an external ultrasonic 
field on the radiation of channeling particles in a crystal was given in 0. In this article 
(see also [Q) spontaneous emission by channeling particles in a silicon crystal bent 
by an acoustic wave of a frequency ~ 10 MHz was briefiy considered. The estimates 
carried out in p, Q were based on the so-called dipole approximation for the undulator 
radiation (see e.g. [^). The latter implies that the strong inequality, p ^ l-K^ajX <^ 1 
is valid, so that all the undulator emission occurs in the fundamental harmonic. Such 
an assumption leads to a considerable narrowing of the parameters of the AIR radiation 
and, also, disregards the possibility to generate the emission in higher harmonics. 
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In papers ^ the spontaneous channeling radiation in the presence of a low- 
amplitude ultrasonic wave was considered by means of classical electrodynamics. In this 
limit the AW amplitude is much less than the interplanar spacing, a d, resulting, thus, 
in low values of the curvatures of the acoustically bent channels. This inequality allowed 
the authors of 0, ^ (see also |^) to consider the phenomenon of resonant enhancement 
of the photon yield due to the coupling of two mechnisms of the photon emission, the 



channeling radiation |2^, and that induced by the AW. 

In 10 both the spontaneous undulator emission and that of a free electron laser 
type was investigated in the case of a relativistic positron beam channeling through a 
periodically strained lattice. In this paper the main attention was paid to the radiation 
emitted in an undulator based on a solid state superlattice (see also p). By using 
the formalism of classical electrodynamics the authors of performed a self-consistent 
treatment involving the wave equation for the electromagnetic radiation field and the 
kinetic equation for the positron distribution function. The important result, which is 
the gain coefficient for the forward radiating field was obtained analytically in a closed 
form. Numerical estimates of the gain in the case of a resonant coupling of the channeling 
radiation and that arising from the superlattice based undulator were presented. 

In there is also a short comment on the possibility of using an acoustically bent 
crystall as an undulator. In our opinion, the analogy between a superlattice and an 
acoustic wave can not be considered as adequate, because the acoustic wave amplitude, 
a, can greatly exceed the interplanar distance d, whereas in the superlattice based 
undulator the inverse condition is assumed, a <^ d. This latter case was analyzed in 

It is stressed in figure |I] and in (|I]) that we consider here the opposite case 
a ^ d. In this limit there is no resonant coupling between the channeling and the 
undulator radiation because the characteristic frequencies of the photons emitted via 
these mechanisms are totally incomparable (see section As a result it is possible 
to distinguish two mechanisms of the radiation formation in an acoustically bent 
crystal and to investigate the properties of the AIR, both spontaneous and stimulated, 
separately from the channeling radiation. 

As we demonstrate below, staying within the regime defined by (|l]), it is possible, 
by tuning the AW amplitude or frequency, to vary significantly the intensity of the 
AIR radiation, to change the patterns of its spectral and angular distributions. These 
characteristics of the radiation, as it is known from the general theory of undulators 
(e.g. [^), are crucially dependent on the magnitude of the undulator parameter p. In 
our paper we provide the analysis of the AIR radiation for both limiting cases, <^ 1 
and p > 1. 

In the case a ^ d it is important to establish the conditions under which 
the acoustically bent crystal may serve as an undulator for the channeled ultra- 
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relativistic particle. A comprehensive analysis of these conditions, and consequently, 
the establishment of the ranges of all parameters in the problem (these are: the AW 
amplitude, frequency and velocity, the type of projectile and its energy, the width of 
the channel and the magnitude of the interplanar field) inside which the AIR process is 
feasible, are studied in our paper. In contrast such a discussion is fully omitted in the 
papers cite above. 

Another point we wish to emphasis is that the quasi-classical approach used in 
our paper is more appropriate than the purely classical treatment utilized in the cited 
papers. In our paper we consider the spectral distribution of the radiation in the whole 
range of the photon energies: huj = . . .e. It is clear that the classical approach is valid 
only in the soft-photon limit, i.e. hu -C e, but is totally inappropriate in the case hu ^ e 
where the quantum corrections (the radiative recoil) must be taken into account. The 
emission of high-energy photons becomes important for energies of projectile positrons 
in the tens GeV range. 

In our paper it is demonstrated, for the first time, that it is feasible to consider, by 
means of the AIR mechanism, stimulated emission within the 10~^ . . . 10° MeV range. 
To establish this interval we carried out a realistic analysis of the influence of the 
dechanneling and the photon attenuation effects on the spectrum of the AIR radiation. 

All the features mentioned above distinguish the results presented below from the 
earlier works 0. 

Finally, let us note that we do not consider the electromagnetic radiation ^ arising 
during the projectile channeling in an ultrasonically excited crystal due to specific 
diffraction of the photons by sets of periodically bent planes. 

The paper is organized as follows: 

In section ^the general formalism of the quasi-classical description of the undulator 
radiation is presented. Section || is devoted to a detailed analysis of the conditions which 
allow to consider the acoustically bent crystal as an undulator. Numerical results for 
the angular and spectral distributions of the undulator radiation as well as the energy 
loss due to the AIR are presented in section ^ for various crystals and projectile energies. 
In section ^ the estimates, both analytical and numerical, for the possibility to generate 
the stimulated emission by means of the acoustically based undulator are presented. In 
section ^ we outline the relevant problems which, to our mind, deserve to be thoroughly 
investigated. 

2. Quasi-classical description of the undulator radiation 

Within the framework of the quasi-classical approach the distribution of the energy 
radiated in given direction by an ultra- relativistic particle (of a spin s = 1/2) and 
summed over the polarizations of the photon and the projectile, is given by the following 
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expression 

dEUn) = (£)'dk I dti dh e-''^(*i'*^) fit,,^) (4) 

Here dk = c~'^uj'^dujdQn, n = ck/cj is the unit vector in the direction of the photon 
emission, c is the hght velocity and q is the projectile charge. The functions ^9(^1,^2) 
and f(ti,t2) are defined as follows 

V9(ti,t2) = ^1 - ^2 - -n- (ri - rs) (5a) 
c 

= i{(l + (H-")')(^-l)+^} (5*) 

The notations used are ri 2 = r(ti 2), Vi 2 = v(ti 2), with r and v standing for projectile's 
radius vector and velocity, respectively. 

Expression looks almost like the classical formula |5T|, although with 
quantum corrections: 

iO > iO = u= — , (d) 

£ — huj e — hu 

which take into account the radiative recoil. 

Let us now obtain a general quasi-classical expression for the radiated energy, 
dEi^{n), in the case of a particle moving in an undulator which contains periods. 
The radius vector and the velocity of the particle in an undulator satisfy the conditions 

r(t + r) = r(t)+Vor, v(t + r)=v(t), (7) 

where T is the time interval during which a projectile passes one period of the undulator, 
and the quantity vq 

vo = ^/%(t)dt (8) 

is the mean velocity along the undulator axis. 
Following and introducing the quantities 

= l^dt^ exp m(t)) , v^T= r exp (i$(t)) , <l>(t) =uj'(t- ^L^li^] (9) 
Jo c Jo \ c J 

expression can be presented in the form: 

where the function D{ri) and its argument are 

smTcr] I ujq \ c J 
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respectively. The parameter ujq = 2tt /T is for the undulator frequency. 

Expression (p!OD together with the definitons ( pT]) clearly exhibit the features typical 
for the undulator radiation |3^. The radiation intensity is proportional to the square 
of the total number of the undulator periods, A''^, reflecting the coherence effect of 
radiation. 

The spectral and angular dependence of the radiation are determined mainly by the 
function D{r]), which is well-known in a classical theory of diffraction [^. This function 



has sharp maxima at the points rj = K = 1,2,3 .. .. It is most clearly seen in the case 
of infinite iVu: 

hm N^Div) = Y.6ir^-K) (12) 

For < oo the function D{ri) has main maxima in the points rj = K = 1,2 .. . where 
its magnitude D{k) = 1. In the interval f] = [K, K + 1] D{t]) has minima in the points 
?7™™ = K + m/N^, m = 1, . . . , — 1 where D(?7™™) = 0. The local maxima are located 
at 7]^^^ {m = 1, . . . , — 1) which can be found from the equation A^u tg nriN^ = ctg tti]. 
The values D(?7™^^) rapidly decrease with m and are much smaller than 1. For m = 1,2 
the corresponding values are D{rif^^) ^ 1/22, Dirif^"^) ^ 1/62. The width A?7 of the 
main maxima of D{ri), estimated as /S.K = 2(77™™ — K), is equal to AK = 2/N^ and 
does not depend on K. 

Since dE^{n) oc D{ri) the frequency- angular distribution of the radiation is 
represented by the sets of harmonics u'j^ (each of the width Au;^ independent on K) 
which are effectively emitted at the angle 9 (measured with respect to the undulator 
axis) 



ujqK a ,/ 2 uj'j^ 



<-77^. ^'-W^ A' = 1.2,3... (13) 



c 



Let indices || and _L indicate, respectively, the parallel and the perpendicular (with 
respect to the undulator axis) components of the projectile radius-vector and velocity: 

r(r) = r||(r) + rx(r), v(r) = V||(r) + vx(r), r = [0,T] (14) 

For 7^1 the following relations are valid with the accuracy up to the 7"^ terms: 




c 

where ^ = T'^ /o^vi(t)dt. 

By inserting (|15D into (|^) and introducing the dimensionless variable %p = ujqt, one 
obtains 

<l>(r) = fUp) =r]'4) + —cose - — n ■ yAt) (16a) 

ujq 2 c 
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{16b) 



and in the following expression for dEi^{n): 



du dflr, 



47r2c ^ 



7^1 + 



u 



2(l + u) 



Relol 



;i7) 



where the quantities I^, Iq and In stand for the integrals 



lo 



2tt 



d^, III 







d^, 



2n 

c 



— e^^^'^M^/^ . (18) 



Within the same level of accuracy as in ([T^) the parameter rj from (|Tl|) reads 




The quantity p, introduced here, is called the undulator parameter and is defined 
by the following relation 



V 



2 



V' = 2Y -T (20) 



c 



The parameter 6*0 is related to p and 7 through 

^« - 5^ - ? (21) 

and has the following physical meaning: in the case ^ 1 (so-called, non-dipole limit), 
60 defines the cone, along the undulator axis, < 6' < ^o, into which the radiation 
is emitted. For < 1 (the dipole case) the emission effectively occurs into the cone 
< ^ < 

The formalism sketched above will be now applied to a planar harmonic undulator, 
in which the quasi-periodic trajectory (|^ lies in a plane. Let the undulator axis coincide 
with the 2;-direction and let its plane be the (?/z)-plane. In this type of an undulator 
the particle trajectory is purely harmonic and is defined as 

y{z) = a sin kz (22) 

We assume that the parameters a and k satisfy the condition ^ = fca ^ 1. In the next 
section, where the undulator based on the acoustic wave transmission in a crystal is 
described in detail, it is demonstrated that this condition is fulfilled. 

For an ultra-relativistic projectile, moving along the trajectory (p2D the following 
formulas are valid for the quantities introduced above (the expressions below are written 
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in the lowest orders in ^ and 7 ^): 



Vo ~ c, 
£■2 



V_L ~ COs{kz) <^ Vii, 



4 



sin 2?/', p 



7]'?/' H — sin 2iIj — — \/2 Oq 6 cos (p sin ■?/' . 

ujQ 4 Wo 



Uq = ck , 



(23 a) 
(236) 
(23 c) 



Here ^ and are the polar and the azimuthal angles of the photon emission, respectively. 

By substituting ( p^3a| - p^3cD into (|T^ and (0) one obtains general quasi-classical 
expressions for the spectral-angular distribution of the radiation in a planar undulator 
?T| . These expressions can be simplified in two limiting cases: <^ 1 and ^ 1. For 



the sake of further reference below the corresponding formulas are exhibited. 



2.1. Non-dipole case: ^ 1 



In the limit N^,p^ ^ 1 the spectral-angular distribution of the undulator radiation is 
given by the following asymptotic expression valid for 6 < 6^. 



d£;(n) 16g2 Nl^f'^ D{r]) 




2 61 



1/3 



cos2/5Ai2(C) 



sin2/?Ai'2(C) 



(24) 



where Ai(^), Ai'(C) are the Airy function and its derivative, respectively, and the 
abbreviations were introduced 



X 



6 



uo p" 

472 ' 

9 cos if 




1 e2gj^2^ 



291 



(3 = 1] i arcsin 6] —3x6 jj^^"^, A 



1 



u 



2(1 + m) 



(25) 



V2 9o 

Other notations used in (p^ have been defined earlier. 

The Airy function and its derivative, both satisfying the conditions Ai((^), |Ai'((^)| ~ 
1, C < 1 Ai(C)) |Ai'(C)| — ^ 0, C ^ 1) define the frequency of the radiated intensity 
maximum, cj^ax- ^' ^ ^max the intensity of radiation exponentially decreases. It 
can be shown that cu^ax ~ Vl"^ ^o- It corresponds to -ft'max ~ P^ which is the largest 
number of the radiated harmonics. 

From (p^)-(|25|) follows (in accordance with general theory of a planar undulator 
pO| , pT| , |36| ) that only odd harmonics {K = 1,3,...) are radiated along the undulator 
axis in the soft-photon (classical) limit, hu/e = 0. Indeed, by putting 9 = and 
M = in (p^)-(p5|) one gets [dE{n)/duj dQn]o=o ^ -0(^7) sin2(7rr7/2). Keeping in mind 
that T] = K ± 1/A^u, one gets i^' = 1, 3, . . .. 
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Equations (^) and ( |25D are simplified considerably ii Tiuj/e = 0. It can 
be demonstrated that in this case (^) represents by itself the classical result for 
d£'(n)/dc(j dfin (see e.g. I^Oj) written in the limit ^ 1. It is worth noting that 
for huj/e = the asymptotic expression (|2^ ) and the exact classical formula written in 
terms of Bessel functions produce close results starting with the value = 2. 

The spectral intensity of the radiation emitted by a projectile moving along the 
path ( P^l ) in an undulator of the total length L can be obtain by integrating (|2^) over 
ip = [0, 27t] and 9 = [0, 6^0], or by deducing it from the general expression given in 
The result is ilTl: 



1 dE g2 cj' / , , 
--^ = -N^p- [Gi{s) + 
L aui c io \ 



1 + 



u 



2(1 + m) 



Go(s) 



(26) 



where the parameter s is defined as s = {u' /uj^'-y'^pf'^^ . The functions Gi,2(s) are equal 
to 

2 



/•oo 

Gi{s) = -2s^/^ dx 



TT — arccos 1 — 



Ai(sa;) 



Go(s) 



a/2 



°° dx 



(chx)^/^ 



Ai'(s(chx) 



2/3-s 



(27) 



In the extreme limit — > 00 expressions (0) and (pTf ) reduce to the analogue 
ones for synchrotron radiation [BT|. 



Dipole case: p^ 

The general quasi-classical expression one finds in |3l|]. We will need here its classical 
limit huj/e = 0. The angular distribution of the radiation is given by the following 
expression 

dE(n) g2 ^222 r /i_(^m2\2>| 

, = o < sin^ + cos^ ^ 7^] } D(t]) 28) 

dcudn^ c(l + (7^)2)2| ^ ^yi + ^^0y)j V ; 

Note, that in the limit p^ <^ 1 the parameter 77 reduces to (see 



V^- {—,+-], 0=[0,1~'] (29) 



Expression ( pH]) is written in the lowest order in p"^ and describes the angular 
distribution of the first harmonic, K = 1. Thus 77 = 1 ± l/N^^ is implied in (^). The 
higher-harmonics emission is strongly supressed in the limit <^ 1 ^ . The spectral 
distribution of the radiation in a dipole case can be obtained by a direct integration of 
{m over if = [0, 27r] and 9 = [0, 7^^]. 
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3. Acoustically bent channel as an undulator 

In this section we formulate the conditions which, in addition to (|lDand (|^), must be 
fulfilled to consider the system "AW + channeled particle" as an undulator. 

Let us consider the case of planar channeling of an ultra-relativistic positively 
charged particle in a channel bent by a monochromatic transverse acoustic wave (AW) 
which is transmitted in the 2;-direction (see figure |l|). The shape of the channel centerline 
is described by the dependence (H). 

Another condition which we impose on the AW is 

^ = A;a = 27r ^ < 1 (30) 
A 

where the quantity k denotes the AW wave number, k = 2tt/X. It will be demonstrated 
below that the condition (^) holds in all cases which are of interest. 



3.1. Motion of the particle in a channel bent by the AW 

Let a particle enter the crystal at the point z = having a small incident angle with 
respect to the (curved) crystallographic plane. The corresponding quantitative criterion 
will be given below. 

In the case of small-angle injection one may use the continuous potential 
approximation to describe the interaction of the particle with the crystallographic plane 
||15| , p5| . This potential, U{p) depends only on the transverse displacement, p, of the 
particle from the mid-plane of the channel ( \p\ < d/2). 

In the case of an ultra-relativistic projectile the classical approach adequately 
describes the projectile motion, so that one may totally disregard the quantum 
corrections to the transverse oscillations of the projectile in the channel due to the 
action of the interplanar potential U. The following estimate proves this statement. Let 
Uo denote the depth of the interplanar potential well. The separation of the neibouring 
energy levels, hojo, of the transverse oscillations is equal, by the order of magnitude, to 
huo ~ /i^g U" /ni'-]' ~ h^qUo/cPm'-y. Hence, the number of these levels, estimated as 



Uo d qUo-f h 

oc Ac = — (31) 



hujo Ac V fnc^ ' i^c 

is large when 7^1, and trajectories can be applied to describe the projectile dynamics. 

The particle's trajectory lies in the (yz)-plane and is governed by the equation of 
motion 

^ = -ggradf/ (32) 
dt 

where p = mvj is the kinematic momentum, 7^^ = (1 — f^/c^ — v^/c^). 
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We assume that the dependence U{p) on p in the case of a bent channel is the same 
as for a hnear one. Below we demonstrate that this assumption is adequate. 

The Hamiltonian of the particle, H = y/p'^c^ + m^c^ + qU{p) does not explicitly 
depend on t, hence, the relation 7 = H/mc^ = 70 (1 + qU/mc'^'-fo) is valid, with 
7o = (1 — f^/c^)~^/^ and v standing for the initial particle velocity. The condition 
qll ^ mc^7o allows to consider 7 as an integral of motion. 

To proceed further and to obtain all the relations analytically we consider the 
harmonic approximation for the interplanar potential, U{p) = Kp'^/2. 

Taking into account (|30|) the transverse coordinate p satisfies the equality p = 
{y — a sin kz) which is valid up to terms ~ ^. With the same accuracy the equation of 
motion (13^) reads 



— {p + z^coskz) = p, -Tw = ^ COS kz (33) 

at m7 dr"^ 7/17 

In lowest order in ^ the system (^) yields z = vt (with v ^ c) and (^ for 
the particle trajectory. Within the linear approximation in ^ the second equation (^) 



defines the deviation from this trajectory through |40 



The term —pv/R{z) represents the centrifugal force due to the curvature of the bent 
channel. R{z) is the curvature radius of the curve (^ and is equal to 



When writing ( |35D we neglected the terms ~ and higher. 

Equation ( P3| ) combined with the inequalities (0) and ( PU] ) justifies the 
approximation made above (see the paragraph after figure ^2]) concerning the equivalence 
of the interplanar potentials in the bent channel and in the linear one. Indeed, the 
deviation of the function U{p) due to the channel curvature is proportional to the ratio 
d/R, which can be neglected: d/R oc {d/X){a/X) <^ <^ 1. 

The general solution of (0) reads 

p{z) = po sm{xz + 0) + — 7^ sin kz, z = vt (36) 

X ~ k 

where we abbreviated = Q'^/ {fri'-fv^). The first term on the right-hand side of 
(p6|) describes the oscillations due to the action of the interplanar force —qU' and is 
responsible for the traditional channeling radiation [Q. The constants po and (p are 
subject to the initial conditions of the projectile injection into the channel. In the case 
when the vector of the initial velocity of the particle is aligned with the tangent of the 
channel centerline one may put po = 0, and, consequently, posin(xz + 0) = 0. 
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The second term in ( p6D corresponds to the oscillations of the projectile in the 
vicinity of the mid-plane due to the action of the (driving) centrifugal force. The 
channeling in a bent channel is stable against the driving oscillations provided their 
amplitude, ak'^ / (x^ — /c^), is less than half the width of the channel: 

< o (37) 



Let us introduce a minimum curvature radius of the channel bent by the transverse AW 

1 / A \^ 1 /V\'^ 1 



^-'"^-^-1,2^; a-fej ^ ^^^^ 

Here V and u are the velocity and the frequency of the AW, respectively. 

Recalling (|l|), the condition (|37D is presented in the form equivalent to (@): 

— < 1 (39) 
R ■ nil' 

Here the notation f/^ax is used to substitute 2n{d/2), and we have put pv ^ pc = e. 

The relation (|39| ) is written in the form which can be used for an arbitrary potential 
U{p). It establishes the condition for a channeling process in a bent channel |T^: the 
maximal centrifugal force, e/R^i^, must be less than the maximal interplanar force 
qU' . 

y ^ max 

The parameter on the left-hand side in (^) defines also a critical incident angle Gc 
25| of the particle injection into acoustically bent crystal: 0c = 0l (1 ~ ^/-Rmin ^^max)' 



where Gl = \j2qUmax/'^l is the critical angle for the linear channel 

Provided (|l|) and (|39D are fulfilled, the projectile, entering the crystal under an angle 
< 0c (with respect to the centerline) will be trapped in the channel and move along 
the trajectory defined by (H). The deviation, (0), of the real trajectory from (|^) can 
be neglected. 

By taking into account the last equality from (|3^) the condition (|39|) may also be 
denoted as 



iy'a<C = ^-' ■ [ — ] ■ ^^^^ (40) 



27r/ Y mc 

and determines the ranges of u = V/X {u is the AW frequency), a and 7 for which the 
channeling process as well as the undulator radiation, can occur for a given crystal and 
crystallographic plane (the parameters t/^iax ^ depend on the choice of the particular 
crystal and the crystallographic plane) and for given projectile type, characterized by 
its rest mass m and charge q. 

Figures ^ illustrate the ranges of u and a in which the channeling process is 
possible for a positron and a proton both of the same energy (as indicated) in C 
(diamond). Si, Ge and W crystals near the (110) crystallographic plane. The solid 
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thick lines represent the boundary v"^ a = C, so that the range of vahdity of (^) 
hes below these lines. To obtain V we used it s relation to the crystal elastic 
constants cn and Cu and the density p: V = ^ (cn — Ci2)/2p. Using the data from 
ones gets the following values of V for C, Si, Ge and W crystals respectively: 
: 11.6 ■ 10^ 4.67 ■ 10^ 2.76 ■ 10^ 2.81 ■ 10^ cm/s. 
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Figure 2. The ranges of i' (in Hz) and a (in cm) in which the channeUing process is 
possible for a positron and proton along the (110) plane in (a) a diamond, (b) a silicon 
crystal, (c) a germanium crystal, (d) a tungsten crystal. The heavy full lines represent 
the boundaries, v'^a = C (see (|^)), for a 5/50/500 GeV projectile as indicated. The 
dashed lines correspond to the constant values of the parameter A^(cm~^) = A^-'^ 
(see also the explanations in the text). The dotted and the chain lines indicate the 
constant values of the undulator parameter p = 2Ti^av /V for a positron and a proton, 
respectively. The values of p, pi,P2 and p^ correspond to the projectile energies 5, 50 
and 500 GeV, respectively. 

The data for qU^^^ (for a positron) are qU^^^ = 12.0, 5.7, 10.4, 43.0 Gev/cm, 
and for the interplanar distance d = 1.26, 1.92,2.00,2.24 10""^cm for C, Si, Ge and W 



crystals, respectively. They were taken from |25, 31|. 

Taking into account the double-log scale in these figures it may be concluded, that 
the channeling condition (^) is fulfilled within wide ranges of the AW amplitudes and 
frequencies, for different types of projectiles and wide ranges of their energies, and for a 
variety of crystals. One can also deduce from figures ^ that within the range of validity 
of (|Dp the condition (|D[) is well fulfilled. 



20 




3.2. The undulator conditions 

Now we formulate the conditions (in addition to the channeling condition established 
above, eqs. (|5I3|)-(PD|)) which allow to consider the system "AW + channeled particle" 
as an undulator. 

Let L denote a crystal thickness. Both the motion of the projectile in the acoustically 
bent channel and the spectrum of the generated radiation are of the undulator type, 
only if 

A < L, (41) 

i.e. if the channeling particle oscillates many times within the length of the crystal. As 
with any other undulator, the suggested undulator is characterized by the frequency. 



cuo, and the undulator parameter, p (see (|20| , p3^) . These quantities depend on the 



AW amplitude and the wave length, and on the relativistic factor, and are expressed as 
follows: 

c a av 

cuq = ck = 271 — , j9 = 74 = 'jka = 27r'j— = 27r7— (42) 
A XV 

For further reference let us note the difference in functional form of p for the AW 

undulator and that for the undulators based on magnetic fields (both helical and planar). 

In the latter case the undulator parameter is equal to (see eg. |^) pb = qB\B/27imc^, 

where B is the amplitude value of the magnetic induction and Ab is the period of the 

magnetic field. The quantity pe is independent on 7 whereas p from (^) is linearly 

proportional to it. 
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In figures ^ the dotted and the chain lines indicate the constant values of the 
undulator parameter p. Note that each dotted and chain hne corresponds to different 
values of p depending on the energy of the projectile. Values of p corresponding to the 
energy 5/50/500 GeV are shown in figure 2 in the vicinity of each line. The broken 
lines correspond the constant values of the parameter both for positron and proton, 
which is defined as the number of the AW periods per 1 cm: N (cm~^) = A~^. Figures 
demonstrate that, within the range of validity of the channeling condition (^), the 
parameters p and N vary over wide ranges: N = 1...100, p = 0.1 — 500 for a positron 
projectile and p = 0.001 — 0.1 for a proton. The upper limiting values of the p's are 
larger by more than an order of magnitude than those accessible in the undulators based 
on the motion of charged particles in periodic magnetic fields ^ . 

The undulator parameter p is proportional to a (see (^21) ). In this connection the 
question arises of the validity of the assumption that a = const during the time of 
flight, T = L/c, oi the particle through the crystal. Hoevwr, for a running or standing 
AW the amplitude is time dependent. To formulate a quantitative criterion let us 
consider the crystal bending by a standing AW. In this case the constant value a must 
be substituted with d(t) = ao cos(27r z/t). For time scales r <^ l/u, one gets the following 
relation for the deviation, Aa = ao — a{t), of a{t) from its amplitude on the timescale 
r: Ad/tto = 27r^(rz/)^. Thus, thie time dependence of the AW is negligible provided 

-It)^©^- (43, 

This condition allows to disregard the variations in the shape (|^) of the channel on the 
timescale r. 

To estimate the corresponding values of L let us consider u < 100 MHz (see figures^. 
Then one gets L ^ 70 cm which is more than well-fulfilled for any realistic L-value. 

Finally, let us demonstrate that the characteristic frequencies of the undulator 
radiation due to the particle motion along the trajectory (^ and those of the ordinary 
channeling radiation are well separated, provided conditions (|I]) and (^D]) are fulfilled. 

The ordinary channeling radiation is due to the transverse oscillations (the first term 
in (^6|)) of the particle inside the channel under the action of the interplanar force. The 
characteristic frequencies of this radiation may be estimated as : Uc ~ 7^x c, where 
X = {qU'^Jedfl' (see (1) and (H). 

The characteristic frequencies of the radiation formed in the acoustically based 
undulator are given by ~ = 7^ci;o (see (^2D). Then, the ratio oj^/uJc is equal to 

J— - « 1, (44) 

provided the conditions (|ip and ( PPD are fulfilled. Indeed, for a ~ lOci and -Rmin ~ 
10 e/qU^^^ one gets Un/uJc ~ 10"^ 
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This condition allows to distinguish two mechanisms for the radiation formed in the 
channeling process in an acoustically bent crystal, and to consider all the radiation in 
the spectral region -C cUc as formed via the undulator mechanism. 

We note here, that the inequality is fulfilled in the whole range of the AW 
parameters consistent with the channeling condition (^), only if d < a (eq. (|I])) is valid 
too. In the opposite case, a < d §|, 0], which can be realized in the undulators based 
on a superlattice or/and on a low-amplitude AW ^ the frequencies and 
ujc can be of the same order of magnitude and, thus, two mechanisms of the radiation 
formation must be treated simultaneously. As mentioned several times already, in the 
present paper only the case d < a is considered. 

When discussing the motion of the particle in a bent channel and the undulator 
conditions, we ignored so far the effects of dechanneling and the photon attenuation. 
The influence of these effects on the undulator radiation will be discussed in the next 
sections. 

Let us summurize our results so far: Provided conditions (|l|), (|4y, ^) are fulfilled, 
the system "ultrarelativistic charged particle + acoustically bent crystal" represents 
a new type of planar undulator, and, consequently, serves as a new source for high- 
frequency undulator radiation of high intensity, monochromaticity and a particular 
pattern of the angular-frequency distribution. The latter as it will be illustrated in 
the next sections. The validity of the classical treatment of the particle motion in a 
bent channel allows to apply the quasi-classical approach, described in the previous 
section, for calculating the characteristics of the undulator radiation. 

4. Numerical results for dE(n)/duj dQ^ and dE/du 

In this section we present numerical results for the spectral and the angular distrbutions 
of the spontaneous undulator radiation formed during a passage of an ultra-relativistic 
particle through an acoustically bent crystal. The calculations are performed by using 
the formalism outlined in section |^. We resctrict ourselves to the case of positron planar 
channeling. 

4.I. The angular distribution of the radiation 

We start with the illustration of the general features of the angular distribution of the 
radiation formed in an acoustically based undulator. 

Figures ^ and ^ correspond to the case of a 50 GeV positron channeling along (110) 
plane in a carbon crystal of thickness L = 1 cm. For chosen values of the AW amplitude 
a = 2.35 ■ 10^^ cm and frequency z/ = 20 MHz, the corresponding magnitudes of the 
total number of the undulator periods, the undulator parameter and the fundamental 
harmonic energy are as follows: = 17, = 10.1, hui = 9.9 MeV. In a diamond 
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crystal the spacing for the (110) planes is d 
well fulfilled. 



1.26-10 ®; therefore the condition 
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Figure 3. The spectral dependence of the radiated energy for a 50 GeV positron 
channehng in a diamond. The polar angle of the emission 9 = Q° . Other parameters 
are as indicated. See also the explanations in the text. 



Figure ^ represent the spectral distribution of the radiation emitted along the 
undulator axis, (di?(n)/da; df2n)6i=o° (^^^ (ID)- Each peak corresponds to the 
emission into the odd harmonics ^ the energies of which are deduced from 
the relation 

= -T^ (45) 

which is obtained from (|l^) and (|15D by putting = 0°. 

It is seen that all harmonics are well separated: the distance between two 
neighbouring peaks is Ihiox ~ 20 MeV whilst the width of each peak is huJi/N^ 1.2 
MeV. In the vicinity of each main maximum the shape of the spectral distribution is 
governed by the behaviour of the function D{ri) defined in (pHI). 

Two 3D plots in figures ^(a) and ^(b) allow to visualize the spatial behaviour of the 
angular distribution of the radiation (at fixed frequency) as a function of the azimuthal, 
ip, and the polar, 6', angles of the photon emission. In these figures the quantity 
Ig (d£'(n)/d(^ci;) dfin) (measured along the 2;-axis) as a function of the azimuthal, (/?, 
and the polar, 9, angles of the photon emission is plotted for two central harmonics. 
Figure ^(a) represents the distribution of the radiation in the odd harmonic Kujk ~ 71.4 



Figure 4. The spatial distribution of tlie radiation emitted in: (a) the odd harmonic 
with K — 7, (b) the even harmonic with K ^ 6. Other parameters are as in figure ^. 
The undulator axis lies along the z direction. The undulator plane is yz. The x 
and y axis are scaled with respect to the dimensionles parameter 9/9o- See also the 
explanations in the text. 
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MeV with K = 7, figure ^(b) corresponds to the emission in the even harmonic 
huK ^ 61.2 MeV, K = 6. 

The undulator axis hes along the z direction, and yz is the undulator plane. The 
positive y direction corresponds to = 0°. In these figures the dimensionless variable 
6/6o is used to characterize the distribution in of the radiation with respect to the polar 
angle. In the considered case the quantity 6*0, defined by (|2l|), is equal to 9o = 23 /irad. 

Figures 0(a,b) illustrate general features intrinsic to the planar undulator radiation 
in the case ^ 1. We first note that the intensity of the radiation in the odd harmonic is 
governed by a powerful maximum in the forward direction, whereas there is no radiation 
in even harmonics for = 0° The latter reaches its maximal values in nearly forward, 
but not in forward, directions. Apart from the intensity peak in the forward (for K 
odd) or nearly forward {K even) direction, the radiation in a particular harmonic is 
emitted in a wide range of the polar angles. In both of the figures there are several 
clearly distinct non-forward peaks in which the radiation intensity reaches the maxima, 
although the magnitudes of dE{n)/d{huj) dfl^ in the maxima at 6* 7^ 0° rapidly decrease 
with the polar angle (recall the log scale along the ^;-axis). 

Another feature to be mentioned is the abscence of the axial symmetry in the angular 
distribution shape. More specifically, the radiation emitted within the undulator plane 
is concentrated in the cone 9i\ ~ = p/(v^7), whereas the cone angle for the photon 
emission in the direction perpendicular to the yz plane equals 9± ~ I/7. Thus, the ratio 
9\\/9_L ~ p/ \/2 characterizes the above mentioned asymmetry in the angular distribution 
with respect to the atimuthal angle ip. This ratio is equal to 2.2 in the considered case. 
This peculiarity most clearly exhibits itself in the shape of the spatial distribution of 
the radiation in the odd harmonics in the forward direction, see figure ^(a). 

By changing the parameters of the acoustically based undulator (which are the 
AW amplitude and frequency, the type of a crystal, and the relativistic factor 7 of the 
projectile positron) it is possible to vary the intensity of the radiation, the energy, width 
and the number of emitted harmonics. To illustrate this statement in figures ^(a,b,c) we 
present the dependencies of the spectral distribution of the radiation emitted at = 0° 
for a 2, 5 and 10 GeV positron channelling along (110) planes in Ge and Si crystals. In 
addition to the parameters indicated in the figures there are other ones, i.e. the crystal 
length, the number of the undulator periods, the undulator parameter, the fundamental 
harmonic energy and the width of the line. They are the following: 

Figure §(a): L = 0.1 cm, = 36, = 3.25, hcui = 53 keV, hAuj = 2.9 keV; 

Figure |(b): L = 1 cm, A^^ = 43, = 6.9, hcui = 0.23 MeV, hAuj = 10.5 keV; 

Figure |(c): L = 1 cm, A^^ = 43, = 27.8, htui = 0.27 MeV, hAuj = 12.7 keV. 
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Figure 5. The dependence of h^^ dE{n)/duj d^ln at 9 = 0° on the photon energy 
Uw calculated for a positron channeling in an acoustically bent crystals. The types of 
the channels, the positron energies and the AW parameters are as indicated. Other 
parameters are enlisted in the text. 
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Figure 5. (continued) 



4-2. The spectral distribution of the radiation 

For large undulator parameters, p > 1, the spectral distribution of the radiation emitted 
by a charged projectile channeling in a crystal bent by an acoustic wave is obtained from 
eqs. ( PBIPTD - Using these formulae, we calculated the spectral distributions for different 
parameters of the AW and for various crystals. The results of the calculations are 



presented in figures 

Figures |§(a,b) exhibit dependencies of dE / d{hu) / L / (in cm~^) for e = 50 and 
e = 500 GeV positrons moving along the (110) channels in a diamond and a tungsten 
crystal. The AW frequency is fixed at z/ = 10 MHz. In these figures each curve 
shows the cu-dependence obtained at a particular value of the AW amplitude and the 
corresponding value of the undulator parameter p, as indicated. The lowest values 
for the AW amplitude are chosen to be consistent with the condition (|l]), where the 
parameter d is equal to 1.26 A for a diamond and 2.24 A for a tungsten (110) channels 
25| . The largest a- values correspond to the relation e/^R^i^eU^^.^^) ^1/6 (see (|39|) ). 

These figures demonstrate that at a certain frequency u the intensity of radiation 
can be varied in a wide range by altering the AW amplitude, which is proportional to 
the value of the undulator parameter as a = p/ (7^), or by choosing the type of a crystal. 
The latter option affects the value of the undulator parameter via its dependence on the 
sound velocity V: p = {2'Kv/V)a. Therefore, for equal a, v and 7 the AIR intensity, 
being proportional to p (see (^6])), is approximately {Vc/Vw) ~ 4 times higher for a 
tungsten crystal as compared with that for a diamond. 

The radiation intensity can also be varied by changing the AW frequency at 
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Channel: C, (110) 

a=1.3 nm, p=l 

• • a=27 nm, p=13 

c a=52 nm, p=26 

• » a=77 nm, p=38 

« a=103 ran, p=51 

Channel: W, (110) 

a-2.2 nm, p=5 

• • a=12 nm, p=25 

c o a=21 nm, p=46 

• — — ♦ a=31 nm, p=66 
c o a=40 nm, p=87 




Figure 6. The spectral intensity (per one period) of tlie AIR, dE / d(huj) / L / Nu, 
emitted by: (a) a 50 GeV positron, (b) a 500 GeV positron, channeling along the (110) 
plane in C and W crystals calculated for the fixed AW frequency v = 10 MHz and for 
various values of the AW amplitude a and the corresponding undulator parameters p 
(as indicated). Further explanations are given in the text. 
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(b) 



e =500 GeV 



V = 10 MHz 



Channel: C, (110) 

a=1.3ran,p=6 

• • a=3.5 ran, p=17 

a o a=5.8 nm, p=28 

« » a=8.0 nm, p=40 

c o a=10.3 nm, p=51 

Channel: W, (110) 

a=2.2 nm, p=49 

•— — • a=2.7 ran, p=59 

o a=3.1 ran, p=68 

♦— — ♦ a=3.5 ran, p=77 
c a=4.0 nm, p=87 




Figure 6. {continued) 
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e =5GeV 
a = 10 nm 
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Channel: Si, (110) 

10.0 MHz, p=1.2 

• • 26.0 MHz, p=3.0 

c = 34.0 MHz, p=3.9 

< » 42 MHz, p=4.8 

Channel: Ge, (110) 

10.0 MHz, p=1.9 

— — 22.8 MHz, p=4.4 

o o 35.6 MHz, p=6.8 

Channel: W, (110) 

10.0 MHz, p=2.2 

23.3 MHz, p=5.1 
'>--o 36.5 MHz, p=8.0 
< ♦49.8 MHz, p=10.9 




Figure 7. The spectral intensity (per one period) of the AIR, dE / d(huj) / L / Nu, 
emitted by: (a) a 5 GeV positron, (b) a 50 GeV positron, channehng along the (110) 
plane in Si, Ge and W crystals calculated for the fixed AW amplitude a = 10 nm and 
for various values of the AW frequency u and the corresponding undulator parameters 
p (as indicated). See also explanations in the text. 




Figure 7. {continued) 
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Figure 8. The spectral intensity (per one period) of the AIR, dE / d{'huo) / L /Nu, 
emitted by a positron of various e (as indicated) channehng along the (110) plane in 
Si and W crystals calculated for the fixed AW amplitude a = 10 nm and frequency 
i/ ~ 10 MHz. The corresponding undulator parameters p are indicated. 



some constant value of a. This possibility is illustrated by figures 0(a,b) where the 
dependencies of dE /ditiui) / L/N^ on for e = 5 and e = 50 GeV positron channeled in 
Si^ Ge and W crystals are presented. 

The spectral distributions calculated at different energies e of a positron channeling 
in Si and W crystals are compared in figure |. 

The features to be discussed in connection with all these figures are the pattern of 
the spectral distribution shape and its transformation with the increase in e. It was 
noted above (see the discussion below (H)) that for u' ^ tu^ax ~ the intensity of 

the undulator radiation (both the angular and the spectral ones) exponentially decreases 
to zero. Within the frame of the quasi-classical approximation the photon frequency 
u is related to the quantity u' through the first equality in (||). Thus, the intensity of 
radiation effectively falls off for the photon energies satisfying the relation 



e 1 + e / p'j'^ hujQ 1 + a/Xcp'^ 

Here Ac = Ti/mc = 2.42 ■ 10~^° cm is the Compton wavelength of the positron. To obtain 
the equality in ( ^61) we made use of p = •yka and uq = kc. 

In the case of (a/Ac)p^^ ^ 1, which is achieved by choosing relatively low values of 
the parameters 7, a or z/, the spectrum of radiation is concentrated in the region of low 
photon energies, hu/e <^ 1 whereas for hu/e < 1 it rapidly drops to zero, as follows 
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from the estimation (|i6|). Such a behaviour most vividly is illustrated by figure |^(a) 
(see also the curves for the C (110) channel in figure |^(b) and for the Si (110) and Ge 
(110) channels in figure 0(b)). 

By increasing 7 the parameter (a/Ac)p~^ becomes < 1, leading to a prominent 
increase in the intensity of radiation in the high-energy photon region, huj/e < 1. 
The curves for the W (110) channel in figure ^(b) and the curve for a 250 GeV 
positron channelling in W (110) (see figure ^ clearly exhibit the transformation of 
the spectral distribution shape: it becomes fiatter, thus refiecting that, except for the 
regions huj/e — > and huj/e — > 1, the probability of the photon emission is almost 
independent on the photon energy. 



4-3. The energy loss due to the AIR 

The stability of spontaneous undulator radiation formed during the passage of the 
projectile through an acoustically bent crystal of total length L to a great extend is 
subject to the condition Ae/e -C 1, where Ae: is the absolute decrease in the projectile 
energy, e, due to various inelastic processes. In particular, the above-written inequality 
implies, that the energy loss, Ae, due to spontaneous AIR only must be small compared 
with e. 

To estimate this quantity in the case ^ 1 it is necessary to integrate dE/duj from 
(p6|) over the interval huj = . . .e. On the other hand, to obtain Ae for < 1 it is 
sufficient to use the classical expression (^) for the spectral and angular distribution 
of the radiation and to carry out the integrations over the angles ip and 6 and over the 
photon energies. 

In figures |(a,b) the dependencies of the relative energy loss, Ae/e, on the energy 
of the projectile e is plotted for a positron channeling in (110) channels for C, Si and 
W crystals and for various AW frequencies as indicated. The curves in figure ^(a) 
correspond to the AW amplitude a = 2 nm, which is equal, by the order of magnitude, 
to lOd for all crystals considered. In figure ^(b) a = 10 nm. For all curves in these 
figures the AW frequency and amplitude are subject to the strong inequality z/^a ^ C 



(see eq. (40)). Both figures correspond to a crystal length L = 1 cm. The total number 
of the undulator periods, is related to L, v and the AW velocity V according to 
N^ = L- Vjv. 

Except for the case < 1, it is hardly feasible to establish by explicit analytical 
means the ranges of parameters inside which the inequality Aeje <C 1 is valid. 



Indeed, expressions ( p6D and (^) exhibit quite complicated dependence of di^/do;, and, 
consequently, that of Ae on £, a and v. In addition to these the energy loss is dependent 
(indirectly, through the conditions (|5), (^) and (PD) on the parameters of the crystal, 
i.e. on ?7^ax ^^"^ ^ ^^"^ length L. 
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positron energy, GeV 



Figure 9. The relative energy loss, /S.e/e, due to the acoustically induced radiation 
versus the energy of a positron channeling along the (110) plane in C, Si and W 
crystals calculated for the fixed AW amplitude: (a) a = 2 nm, (b) a = 10 nm, and 
for various values of the AW frequency u (as indicated). The curves correspond to the 
crystal length L = 1 cm. 
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Figure 9. {continued) 
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Thus, for any particular set of the above-mentioned parameters, the magnitude of 
the energy loss must be calculated numerically. Nevertheless, figures ^ allow to estimate 
the range of validity of Ae/e ^1. It is seen that for a = 10 . . . 100 d and L < 1 cm 
there are wide ranges of positron energy [e from < 1 GeV up to several tens of GeV) 
and AW frequency {u from several MHz up to several tens of MHz) for which all the 
undulator conditions ^3|) are fulfilled and, simultaneously, the relative energy 

loss due to the acoustically induced radiation is negligibly small. 



5. Stimulated emission 



In this section we discuss the possibility to generate stimulated emission of high energy 
photons by means of a bunch of ultra-relativistic particles of length Lb moving in a 
channel bent by a transverse acoustic wave. 



The mechanism of the emission stimulation is illustrated in figure |T0|. The photons 
emitted in the nearly forward direction at some maximum/minimum point of the 
trajectory by a group of particles of the bunch stimulate the emission of the same 
photons by another (succeding) group of particles of the same bunch when it reaches 
next maximum/minimum. 



emitted photons 



3-:- 



the bunch of parti 

Figure 10. Mechanism of the radiation stimulation in the AW based undulator. 
Photons, which are emitted in the forward direction by the positron bunch (dark 
block) as it passes the undulator maximum (marked with A), stimulate the emission 
by the same bunch in the vicinity of the next maximum (marked with B). It is assumed 
that the length of the bunch, Lf, , its velocity, v, and the undulator period A, satisfy the 
condition (1 — v/c)X <C Lb, which means that the photon slippage against the positron 
bunch during one undulator period is much less than the latter. 



This scheme implies that the stimulation is due to the motion of the same bunch 
along the trajectory (^). Therefore, the bunch velocity v, its length Lb and the length 
of the crystal L must satisfy the condition 

L L- Lb V Lb , - 

— > , or - > 1 (47) 

c V c L 

which means that the electromagnetic wave does not slip against the bunch of particles 
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by more than Lb on the lenghtscale L ^ L\y and, thus, the amphfication occurs during 
each of the undulator periods of total number = L/\. 

In the theory of free-electron lasers (FEL) this principle is called "Self-Amplified 
Spontaneous Emission (SASE)" P8| , ^ and usually is referred to the FEL operation in 
the high gain (or collective instability) regime |Q, ^ 



An accurate treatment of the evolution of the photon pulse in the SASE mode implies 
that the space charge fields, energy spread of the particles in the beam and diffraction 
effects are taken into account In the case when the beam undergoes channeling 

the picture must be additionally complicated by considering the dechanneling dynamics 
as well as the photon attenuation jS^ . 



We do not aim to fulfil this program in the present paper. Instead, for the suggested 
undulator, we shall present an estimate of the gain exploiting the general relationship 



between gain and spontaneous spectrum of undulator radiation |3^, |55[. Thus, we 
disregard the dynamics of the particles in the undulator and consider the case of a 
monoenergetic beam of independent particles. Additionally we assume that there is 
no beam divergence. The effects of dechanneling and photon attenuation are discussed 
qualitatively on a phenomenological level. 

5.1. General expression for the gain 

It is our aim to establish the range of parameters (which are: the AW amplitude and 
frequency (wavelength), the type of the particle and its energy, and the type of a crystal) 
within which there is a principle possibility to consider the SASE radiation by means of 
the acoustically bent crystal. 

It had been established long ago that in a free-electron device the stimulation 



of emission occurs if only the photon frequency is equal to that of the spontaneous 
undulator radiation, ujk (see (0)). 

Let us estimate the gain factor, qk, which defines the increase (i.e. the case qk > 0) 
or the decrease (if qk < 0) per 1 cm in the total number, Nk, of the emitted photons 
at a frequency ujk 

dNK = gKNKdz, (48) 

due to stimulated emission (absorption) by the particles in the beam. We shall do this 
within the frame- work of low-gain approximation (see eg. [Q). Also we restrict our 
consideration to the case of the stimulated emission of soft photons, huj <^ e, and, hence, 
u' = u (see (ID). 

The general expression for the quantity qk is 

gK = n [ae{€, e - hux) - (y^{e, e + hux)] (49) 

Here cre(£,£ — hux) and o"a(£,e + hux) are the cross sections of, correspondingly, the 
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spontaneous emission and absorption of the cui^-photon by a particle of the beam, n 
stands for the volume density of the beam particles. By using the known relations 
between the cross section of the photon emission/absorption and the spectral- angular 
intensity of the emitted radiation and taking into account the relation huj <^ e, one 
gets the following expression for the gain, corresponding to the increase in the number 
of photons of the frequency within the interval uk ± Au; emitted in the cone AQk (see 
below) with the axis along the the z-direction: 



9k 



-(2vr) 



n ■ 



K 



de 



dE 



from 



The width Acu one gets 
(2/iV,) (47W(2+P')). 

Making use of ([13|) (see also ( [T9| 

277(1 + 2p-2) 



dujK d^l 
(H by 



(50) 



substituting uj with u, 
one gets the estimation of the quantity AQk'- 



AQ 



K 



K 



(51) 



It defines the maximal solid angle along the undulator axis inside which the Kih 
harmonic is well resolved. The quantity 9q is given in (pT]). Expression (51) is valid 



for both ^ 1 and <^ 1. In the latter case one may set (1 + 2p~'^)6q 
For the energy emitted in the forward direction [6 = 0) the formulas 
yield: 

,2 



7" 



and 



dE 



dui dQ 

dE 
duK dfl 



c 



for p <^ 1 



c p 



for ^ 1 



(52 a) 
(526) 



Expression ( |52 o|) is denoted for the fundamental harmonic, Ui = 2c(Jo7^ in accordance 
with (pSl). In ( |52£| ) we assumed that K <^ -K"max ~ 



and took into account that in this 



case the argument ( (see (0)) of the derivative of the Airy function from (p^ is small, 
so that 16Ai'2(0) ^ 1 pf. 

When carrying out the derivative of [dE / dui dfl\0=o with respect to e the main 
contribution comes from the term dD (77) /de = {d'r]/de)-{dD{r])/d'r]). To obtain maximal 
positive gain it is necessary to choose that slope of D{ri) where it is negative. Then the 
extremal value of this derivative is equal to 



dPjr]) 

de 
dD{7]) 

de 



drj/de 
drj/de 



dD{ri) 

di] 
dD{7]) 

dr] 



— ■ — tor p <^ 1 

e 2 

— for p ^1 



ep^ 



(53a) 
(536) 



The first factors on the right-hand sides of ( |53a| ) and ( p3fc| ) correspond to (d?7/de) and, 
as it is seen, their magnitude depends on p"^. The last factor equals to the 



36 



maximum value of the derivative {dD{ri)/dr]) with t] lying within the interval K±l/N^. 
We concentrate our attention on these details in order to point out the nature of the 
difference in the p-dependence of gx for the AW undulator and that for the undulator 
based on magnetic field. As it was noted above (see the paragraph after (^21) ) in the 
latter case the undulator parameter is independent of 7 resulting in d?7/d£ = —2K/e for 
both p| -C 1 and p| 3> 1 cases {K = 1 if p| -C 1). In turn it leads to the proportionality 
of the gain to p| regardless to the magnitude of the latter ET^. In our case the gain 
gx, which is equal to 



g2 J\j2 ^ 

gx ( cm"^ ) (2v:f — —V < 1 (54a) 



9K{cm ^ ' ^ rci n for / > 1^ 54^ 

is independent on p if ^ 1. In the above equations Ac = h/meC = 3.9 • 10^^^ cm is 
the Compton wavelength and = e^/meC^ = 2.8 ■ 10~^^ cm is the classical radius of 
the electron, q and m are the charge and the mass of the projectile measured in units 
of the elementary charge e > and electron mass mg, respectively. The quantity A is 
measured in cm and n in cm~^. 

Note the strong inverse dependence of gx on 7 which is due to the radiative recoil. 
The gain is proportional to the factor A'"^, which reflects the coherence of radiation. The 
proportionality of the gain to A means that the increase in A leads to an enchancement 
of the radiation intensity in the forward direction. 

Expressions ( ^4aD and (545) were obtained by using the asymptotic formulae (^) 
and ( p8l) valid for 1 and <^ 1, respectively. Extrapolating ( ^4aD and (|54^ ) to 

the region ~ 1 and considering in (|54&|) the fundamental harmonic, loi ~ A'j^ujq/p'^, 
one recognizes that in the point p = 1 (54&) is within a factor 2~^/^ = 0.79 of (|54q|) . 
Such a discrepancy may be ignored having in mind the approximate character of the 
initial formulae (0) and (pSl) . 

Apart from the absence of the p'^ dependence in case p^ <^ 1 (due to the reasons 
mentioned above) expressions ( ^4q|) and ( |546|) have the same functional form as those 



obtained for the undulator based on the action of periodic magnetic field |E6|, |58|, |57|, |47[]. 



5.2. Estimations of gx for the AW based undulator 

Expressions ( |54q|) and ( |546|) allow to make quantitative estimates concerning the 
possibility of stimulated emission of the high-energy photons. 

Before carrying out this analysis let us note two main qualitative features which 
distinguish the AW based undulator and that based on the action of the periodic 
magnetic field. In the latter case both the particle and the photon beams propagate 
in vacuum. Therefore, there are no strong mechanisms leading to the decrease in the 
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beam density. 

In the scheme considered here both beams propagate in a crystaL Thus, an 
adequate description of the stimulated emission with necessity must account for perticle 
dechannehng and for photon attenuation. 

Random scattering of the channehng particle by the electrons and nuclei of the 
crystal leads to the gradual increase of the particle energy associated with the transverse 
oscillations in the channel. As a result, the transverse energy at some distance, L^, from 
the entrance point exceeds the depth of the interplanar potential well, and the particle 



leaves the channel. The quantity is called the dechanneling length [pq] . 

The dynamics of the photon flux propagation, apart from the stimulation effect, is 
strongly influenced by a variety of processes occuring in a crystal. These are the atomic 
and the nuclear photoeffects, the coherent and the incoherent scattering with electrons 
and nuclei, the electron-positron pair production (in the case of high energy photons). 
All these processes lead to the decrease in the density of the photon flux as it propagates 
through the crystal. 

The rigorous treatment of the stimulated emission in the system "channelng beam + 
photon flux" , with both the dechanneling effect and the photon attenuation accounted 
for, must include solving the coupled equations of the Fokker-Planck form for the particle 
beam |3^ and for the photon beam |6y]. This will be done at a later time. In 



this paper we present the quantitative analysis of the stimulated emission by a beam 
channeling in the AW bent crystal by using phenomenological arguments when treating 
the dechanneling and the photon attenuation. 

Let us first estimate the role of the dechanneling effect. 
To simplify the expressions, we will only consider the stimulated emission of the 
fundamental harmonic 

2 +p2 

and omit the index "1" in the notations Ui and gi. Also, we restrict the treatment to 
the case of a positron channeling, M = Z = 1. 

By expressing the total number of the undulator periods A^u through the undulator 
length, Na = L/\, and using (|54q|) and ( ^4^) , one gets the following estimate for g: 

g = g{L)^Q{p){27if r,,n^ (56) 

where the quantity 6(j9) is defined as 

1 for >l 
for <1 



0(P) = <I 2 . (57) 



In ( 1^ we have omitted the factor 2 = 0.79 ~ 1 in the case > 1. In (^) the 
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argument L in the notation g{L) stresses the dependence of the gain factor on the crystal 
length. 

Provided the dechanneling is neglected, one may unrestrictedly increase g by 
considering larger L-values. In reality, the volume density of the channeling particles 
decreases with the penetration distance, z, and, roughly, satisfies the exponential decay 



law 25 



n{z) = no exp {~z/ Ld{^, R)) (58) 

where no is the volume density at the entrance, -^^(7, R) is the dechanneling length. For 
given crystal and channel Ld{'y, R) depends on a positron energy (relativistic factor) 
and on the curvature radius R. 

As argued in p5[, the decrease of n with z acquires the exponential form (|58|) at a 
considerable penetration depth, z > L^ipf, R). If the crystal length is less than -^^(7, R) 
then the dependence n{z) is strongly influenced by the initial conditions (the incident 
angle, the beam divergency) of the beam impacting on the crystal. Therefore, for 
estimation purposes, we set the length of the crystal equal to the dechanneling length, 
L = Ld{'~f, R), and assume that the beam density for z < L is equal to its initial value 
no. 

For a bent crystal with a constant curvature radius R the dechanneling length 



Lrf(7, R) satisfies the relation ||25 



L,(7, R) = il- x2)%(7, 00) ^ (1 - §)'^d(7, 00) (59) 

where -^^(7, 00) is the dechanneling length of a positron of the same energy in a straight 
channel [R = 00) and 

Rc = ^ (60) 

is the critical (minimal) radius consistent with the channeling condition in a bent crystal, 
"the centrifugal force < the interplanar force" (0). 

In an acoustically bent crystal the curvature 1/R{z) is not constant (see (PS|)). 
Therefore it is natural to consider the mean curvature 1/R which is obtained by 
averaging l/|i?(z)| over the undulator period 



^ = — / k a\ smkz\dz = (61) 

R X Jo TT -Rmin 

with -Rmin from (|38|). 

To proceed further we make use of the fact, that for a given crystal and 
crystallographic plane the reduced dechanneling length, a = Ldipf, c>o)/7) depends on 7 
weakly. Its explicit expression, calculated by using the Lindhard approximation for the 
potential of a planar channel, reads 

256 flTF d 
"^^^ = 9^ 77 In (2.//)-! 
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Here otf = 0.8853Z-i/=^ao and / = 16Z^-^ eV are the Thomas- Fermi atomic radius and 
ionization potential, respectively. Zc is the atomic number of the crystal atoms, and 
is the Bohr radius. The dependences 0(7) for various planar channels are illustrated in 
Figure O. For the LiH crystal we used the average atomic number Zc = 2, and the 



interplanar spacings d for its (100), (110) and (111) channels were deduced from [J3 




Figure 11. Dependence of a straight channel reduced dechannehng length, defined 
as a = Ld{j, 00)/ J (see also (|6^)), on the relativistic factor 7 of a positron channelled 
in (100), (110) and (111) planar channels of C, Si, W and LiH crystals as indicated. 
The full-circle curves correspond to (111) channels, the open-circle curves to (110), the 
curves without circles stand for (100) channels. The analogue curves for the Ge crystal 
(not presented in the figure) almost coincide with those for the C. 



By introducing (|59D and (0) in (^61), and expressing k through x, Rc and a as 
k = x^jTT/2Rca, one gets 

giL,) = (27r)2 y|e(p) r„ n ^ x ■ {1 - x'f (63) 

This expression conveniently exposes the dependences of the gain factor on the projectile 
energy, on the parameters of a crystal, t^^ax ^^^1 a, and on the AW amplitude a. The 
factor X ■ (1 — x'^Y contains the dependence of g{Ld) on the average curvature radius of 
the channel since is defined as = Rd R- 

By analysing ( |63D one obtains the ranges of parameters inside which the process 
of the emission stimulation occurs most efficiently. The factor x ■ (1 — x'^Y reaches its 



40 



maximum value of 0.208 at 

defining, thus, for given 7 and crystal, the optimal value of a minimum curvature radius 
of the channel bent by the AW (recall 



KL = - ^ (646) 

Our analysis of an undulator based on an acoustically bent channel assumes the 
validity of the condition (|I]). It follows then, that the range of the a- values, being 
subject to ( |64q|) ), is restricted from below to some a^i^ ^ d. 

Another condition to be fulfilled follows from (^) and implies that the total number 
of the undulator periods on the lengthscale ^^(7, R) = (1 — x^yLdH, 00) = 

L J x=xo 

0.79 ■ Lrf(7,oo) is large. From (|59|), ( pO| ) and ( |64a| ) it follows then 
0.79xo 07 , 



This inequality defines the upper limit of the AW amplitude: 



The inequahties a^i^ ^ and (|66|) must be satisfied simultaneously, leading to a 
natural condition amin ^ c^max- The latter results in the relation 



, 2 /OJ9xo \'mc2 
«7 >- ^Ar_a_, (67) 



max 



giving the lowest value 7min (which is obtained from (^)) with the equality sign). It 
allows to obtain the non-zero mesh defined by ( |64o|) . 

In table |1| the 7min-values are presented for a positron channeling near various 
crystallographic planes in C, Si, Ge, W and LiH crystals. The data correspond to 
Omin = lOc?, Nmir, = 10. Thc d valucs for C, Si, Ge, W were taken from [B5|, BTl, the 



LiH data were adopted from The values of et/^iax were calculated by using the 

Moliere approximation |]T^ for an interplanar potential at the temperature T = 150 K. 

Thus, for given 7 > 7min the interval [amin, Omax] together with the values of 
u following from ( |64fcD and (^), represent the ranges of the AW amplitudes and 
frequencies whithin which the emission stimulation occurs with maximal efficiency and, 
simultaneously, the physical conditions (|I]), (|39|) and ( PD are fulfilled. For these a and 
u values the expression for g is given by 



[^?(^<.)L=.o = 10.3r,in^-^ (68) 



Here the coefficient 10.3 = (2 • xo(l — x\ 



2\2 

oJ ■ 



41 



Table 1. The values of 7min (see (|6^)) calculated for a positron channelling in (100), 
(110) and (111) planar channels in C, Si, Ge, W and LiH crystals. The minimum 
AW amplitude, amin, is set to lOd (d is the interplanar distance). eJ/max stands for a 
maximum gradient of the interplanar potential. Ld^jmin, R) and R are, respectively, 
the dechanneling length and the average bending radius, both corresponding to 7niin- 



Channel 


d 


eU' 


Tmiii 


R 


Ld{lmin,R) 




A 


GeV/cm 




mm 


mm 


C 












(100) 


0.89 


4.8 


1570 


1.68 


0.184 


(110) 


1.26 


7.4 


630 


0.43 


0.111 


(111) 


1.54 


9.5 


380 


0.20 


0.084 


Si 












(100) 


1.36 


4.8 


1650 
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Let us now define the amplitude G [amin, c^max], for which the gain factor (^) 
reaches its maximum. This quantity strongly depends on the factor Q{p)/^/a, with 
e(p) from (13): 

, . ( for > 1 

^= vrxS (69) 
^ 1^ • for p2 < 1 ^ ^ 

I 2Rr 



The latter relation is obtained by replacing k in p = 'jka with k = x^J n /2Rca (see 



(H^) together with (|9|)). 

It follows from (|69[) that the maximum point is at oq = Omin if > 1 for all a from 
the interval [amin, Omax], and = Omax if in the whole interval of a the inequality <1 
is valid. A third option appears if the curves representing the dependencies -^ka = 1 
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and k'^a 



TT 



/2 Xq i?^ ^ cross in the point 
2Rc 



2 9' 

which hes within the interval \a 
Therefore, the quantity 



(70) 



mm; "^max 



. In this case oq = 7 



ao 



'^min 

a 



if a < amin 

if '^min ^ 0- ^ C^max 
if flmax < a 



(71) 



defines the AW amphtude at which the gain achieves its maximum value which is: 



a 



4.3 Vein — X < 



if a < Omin 

if flmin < a < fln 
if Orviax < d 



(72) 



The first harmonic frequency corresponding to the maximum gain (^), is obtained 
from relation (|55D, by expressing uo and through oq and d. The result for the photon 
energy, measured in MeV reads 



hu (MeV) = 3.3 ■ 10 



-11 



(73) 



Figure |T2| illustrates the dependence of the optimal AW amplitude oq (the diamonds- 
embedded curves) and of hu (the circles-embedded curves) on 7 > 7min in the case of 
a positron channeling near the (100), (110) and (111) planes in a diamond taken at 
temperature T = 150 K. The curves corresponds to the particular values of Omin = lOd 
(see table P and A^min = 10. For other crystals the dependencies are basically the same, 
differing (not radically) only in ranges of the parameters. 

The curves in figure |12| clearly demonstrate, that it is meaningful to consider the 
stimulated emission of high energy photons, hu ~ 0.01 ... 10 MeV, during the channeling 
of £ ~ 0.1 ■ • ■ 10 GeV positrons in a crystal bent by a high-amplitude (a ^ d) transverse 
acoustic wave propagating along the channel axis. 

The above consideration focused on establishing the role of the dechanneling effect 
in the process of stimulated emission. The main restriction due to the dechanneling 
concerns the length of an undulator, which must be less than the dechanneling length, 
L < Ld. This, in turn, results in narrowing the ranges of 7 (eq. (|67D), a (eq. (|66D) and 
v for which stimulated emission can effectively occur. 

As it will be demonstrated below, small values of the undulator length, and, 
additionally, the attenuation of the photon flux in a crystal, both require high positron 




2 3 4 

10 10 10 

7 

Figure 12. The magnitudes of tlie maximum gain acoustic wave amplitude (in A), ao 
(see ([TI])), and the first harmonic photon energy (in MeV), huj (see (|7^)), versus the 
relativistic factor of a positron channeling in diamond (100), (110) and (111) planar 
channels (as indicated). 



beam densities necessary to achieve a noticeable level of emission amplification. From 
the theory of FEL it is known that a model which considers a positron bunch as a 
group of independent particles (which has been utilized above) is adequate provided the 
following inequality is valid |^ P7| : 

tOpT<^^/^ (74) 

where cUp = (47rne^/m)^/^ is the beam plasma frequency, and r = L/c is the time of 
flight of the bunch through the crystal. The physical meaning of (|7^ ) is that r is much 
smaller then the timescale needed to develop collective instabilities within the bunch 
due to its interaction with the undulator field and that of the radiation |^|. Therefore, 
all effects due to the space-charge dynamics during the channeling can be neglected if 
the volume density of the beam particles satisfies the condition 

Recalling that Ld{;~i,R) = (1 — a;o)^7a, and inserting the values for the reduced 
dechanneling length a from figures 5.2, one obtains from (|7SD an estimation for n in the 
range of 7 10^ . . . 10"^: n < 10^^ . . . 10^^ cm"^. 

Next we consider photon attenuation. Its influence can be accounted for in a rather 
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simple way by suggesting that the increase per 1 cm in the total number of the emitted 
photons is given by 



dA^= {g- ^)Ndz 



(76) 

rather than by Here the quantity /i = ^{uj) stands for the mass attenuation 

coefficient. These mass attenuation coefficients are tabulated for all elements and for a 
wide range of photon frequencies. For ttui > 1 KeV and up to 100 MeV the /i-values can 
be found in ref. [|^. Figure |l^ shows the ^{uj) dependences for various crystals. The 
step-like behaviour of the curves for Si, Ge and W is due to the discontinuity of the 
atomic photoeffect cross section in the vicinity of inner-shell thresholds. For photons of 
energy less than 1 KeV the main contribution to the photon attenuation arises from the 



atomic photoeffect process, and the n{u) dependence is easily deduced from ref. . 
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Figure 13. The mass attenuation coefficient, fj,, in dependence on photon energy for 
various crystals, as indicated. The data are taken from [p4|. 



To provide an increase in the radiation intensity the gain must satisfy the condition 
9 > /u(^), (77) 

which leads to further narrowing of the interval of photon energies. The large magnitude 
of for hu < 10 KeV (see figure |13[) makes it almost impossible to obtain stimulated 
emission in this ti;-region even for crystals with small average atomic number. 

To draw a conclusion on what type of a crystal and a plane is more suitable from the 
viewpoint of its efficiency to amplify the spontaneous emission formed in an acoustically 
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based undulator let us make qualitative estimates of the density of a positron beam 
needed to achieve the magnitude of 1 for the total gain. 

The total gain, G, which represents by itself the total increase in the number of 
photons due to the stimulated emission on the lengthscale L = Ld{i,R) , equals 

G={g{L)-ii{uj))L (78) 



Such a definition of the total gain is valid in the low-gain limit, G <^ 1 (see e/g. |3T1 , 
Nevertheless, for estimation purposes, one may use ([78|) in the region G < 1. 

Putting G = 1 and representing the quantity g{L) in the form g{L) = g{L) ■ n (the 
explicit expression for g{L) is clear from (|68|) ), we reach the following estimate for the 
volume density of beam particles: 

n = i±fM^ (79) 

The dependences n versus uj for various energies of a positron beam and for various 
types of planar channels (as indicated) are presented in figures |14[ For each channel 
and for e > ^min = mc^7min (the definition of 7min is given by (|67D and its particular 
values are presented in table |l]) the quantity g{L) can be calculated from (|68|) ), and the 
corresponding uj range was calculated by using the relation 

-7^ 1 

/icj (MeV) = 3.3- 10""^= -— -rr (80) 

where a is defined in (|7Dp. In (^) it is assumed that the AW amplitude a varies within 
the interval a = [amin, Omax], and both a and the AW frequency u are subject to ( |64fc| ). 
The minimum value of the AW amplitude was chosen as amin = 10 d, the quantity 
flmax was calculated from (^) with the munimum number of the undulator periods 
A^min = 10. The magnitudes of the mass attenuation coefficients, fi{uj), were obtained 



by interpolating the data from [54 



For each curve in figures 14 the minimum value of the beam density is achieved at 
a = (see (^l])) with the corresponding magnitude of the photon energy given by ([73|) . 
The only exception is the n{uj) dependence for a 250 MeV positron beam channeling 
along the (100) plane in W (see figure 0(a)). In this case, due to the irregularity in the 
behaviour of fj^{uj) in the vicinity of the Is atomic threshold {lis = 69.5 keV for W), - 
see figure |1^, - the minimum of the beam density is achieved at hu ~ 70 keV. 

Figure 15 illustrates the extent to which the photon attenuation infiuences the 
minimum value of the beam density. For each crystal (as indicated) the dependence 
n{u) from ([79| ) is represented by the curves without full circles. The curves embedded 
with full circles correspond to the function h^uj) = ['^(i^)]^(t^)=o; giving thus the minimum 
values of the volume density which is needed to achieve the total gain G = 1 in the case 
when the photon attenuation is totally disregarded. It is seen that for hu > 1 MeV 
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Figure 14. The volume density needed to achieve the total gain G = 1 for a positron 
beam (for various energies e as indicated) channeling along (a) (100) plane, (b) (110) 
plane, (c) (111) plane, in various crystals as indicated. The curves correspond to (|79|). 
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Figure 14. (continued) 
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Figure 14. (continued) 



the role of the attenuation is neghgibly small for all crystals, whereas for less energetic 
photons it leads to a prominent increase in n. 

Figures |l^ and |15| show that to achieve a noticeable level of the amplification of the 
spontaneous radiation formed in the acoustically based undulator large magnitudes of 
a positron beam volume density, n ~ 10^^ — 10^^, are needed. This is mainly due to 
the dechanneling effect which imposes strong restrictions on the undulator length. 

These values of n, being high enough, are, nevertheless comparable with those 



planned to be achieved within the TESLA project |62 



Therefore, we believe that the scheme described above may be considered as a new 
feasible source for stimulated emission within the photon energies range from tens of 
KeV up to MeV by means of the positron beam channeling in an acoustically bent 
crystal. 



6. Concluding remarks 

Our investigation shows that a crystal, which is periodically bent by a transverse acoustic 
wave, can be used for the construction of an undulator for a beam of ultra-relativistic 
particles channeled in the lattice. The parameters of this undulator can be tuned by 
varying the AW amplitude and frequency, the energy of the projectile, and by using 
different types of crystals and its channels. The suggested undulator can be used for 
the generation of the radiation of high energy photons. 
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photon energy (MeV) 



Figure 15. The volume density needed to achieve the total gain G = 1 for a 500 
MeV positron beam channeling near the (111) planar channels in various crystals as 
indicated. The curves without full circles corresponds to ([79|). The curves with full 
circles represent the ( |79| ) dependences with fi{uj) = (i.e. no photon attenuation taken 
into account). 
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Figure 14. (continued) 
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Also, it is shown that it is meaningful to discuss the possibility to create a powerful 
source of a free-electron laser type stimulated radiation in the energy range of tens of 
keV up to the MeV region. 

We consider our present research as an initial milestone for more advanced 
theoretical treatment of the problem. The goal of this further investigation is to achieve 
more accurate quantitative description of the undulator radiation and the corresponding 
laser effect. In out opinion, on this way the following phenomena must be thoroughly 
considered. 

1. The AIR mechanism for a wide beam. When the wide beam enters the crystal it 
is spread over many channels. Each channel is bent by the AW and, hence, one 
must (when calculating spontaneous and stimulated emission) take into account 
the photon fluxes formed by these sub-beams. In this case the coherence effect in 
the photon emission can be enhanced due to the interaction of the photons with 
the particles moving in many neighbouring channels. In the present paper we have 
investigated the undulator and the laser effects in a single channel. 

2. Kinetics in the AIR problem. The more advanced description of the AIR problem 
must take into account the kinetic behaviour of the particles in the beam as well as 
kinetics of the photon flux. Such a description will treat properly the dechanneling 
mechanism, the attenuation phenomenon and the interaction of photons and 
particles beams. 

3. Dechanneling in a periodically bent lattice. It is necessary to study dechanneling in a 
periodically bent channel in comparison with this effect in a straight channel and/or 
in a channel with constant (or/and slowly varied) curvature. So far nobody studied 
dechanneling in a periodically bent channel (the AW bending). The dechanneling 
length might increase as compared with the linear channel case. The dechanneling 
length plays the crucial role for the laser effect which we discuss in our paper (see 
Section Indeed, increasing of the dechanneling length by factor of 10 results in 
decreasing of the channeling beam density, necessary for achievement of large gain 
G, by the factor of 1000 (see eqs. (56) and (79)). 

Let us stress here that special attention must be paid the problem of dechanneling 
of a wide beam in a periodically bent crystal. Broadening of the channeled beam 
in such a medium, which was not considered before, is important in our problem. 

4. The bunching effect. The effect of the stimulation of photon emission can be 
enhanced by increasing density n of the channeling beam. From the theory of 
free electron lasers (e.g. |6^) it is known that for high n large gain factors can 
be achieved in the regime of collective instabilities in the beam (the bunching 
effect). In our future reaserch we plan to apply the results of the FEL theory to 
the description of the laser-type radiation based on the AIR phenomenon. 

5. Acoustic waves of various configurations. In the present paper we investigated 
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spontaneous and stimulated AIR in an acoustically bent crystal considering 
the most simple example of the acoustic wave, the monochromatic plane 
wave. However, other cases of AW (longitudinal waves, spherical waves, non- 
monochromatic waves and various combinations thereof), interacting with the beam 
of the channeling particles, are worthy to study. By applying the more complex 
acoustic waves for the crystal bending one may construct an undulator with variable 
parameters for the generation of high energy photons in a wide range. The use of the 
crystal-based undulator with variable parameters may result in noticeable increase 
of the gain factor for the stimulated AIR analogously to how it occues in the free 
electron lasers based on the tampered magnetic wigglers @]. 
6. Combination of the ordinary channeling radiation and the AIR. The present 
consideration was focused on the case of the high amplitude AW (see (1)). This 
allowed us to disregard the contribution of the ordinary channeling radiation to the 
total spectrum of the emitted photons. More general approach must include both 
radiative mechanisms, as well as their interference, simultaneously. 

The rigorous treatment of the above outlined problems will be the subject of further 
publications. 
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